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Abstract 

Quantum formalism of Fraunhofer diffraction is obtained. The state of 
the diffraction optical field is connected with the state of the incident op- 
tical field by a diffraction factor. Based on this formalism, correlations of 
the diffraction modes are calculated with different kinds of incident opti- 
cal fields. fnfluence of correlations of the incident modes on the diffraction 
pattern is analyzed and an explanation of the "ghost" diffraction is pro- 
posed. 
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1 Introduction 



Correlations of states play an important role in quantum cryptography t 1 ' 2 !, 
teleportation I 3,4 !, and computation I 5 " 11 ! theory. Correlated states are gener- 
ated usually by nonlinear optical processes^ 12 ! or by the beam splitterl 13 ! . In 
this paper, we consider correlations in interference and diffraction. On the one 
hand, the diffraction or interference modes have some interesting correlation 
properties. On the other hand, correlations of the incident modes has a no- 
table influence on the interference or diffraction pattern, in particular, it is the 
key to the explanation of the "ghost" diffraction^ 14 ] , an interesting quantum 
effect. Interference can be regarded as a special case of diffraction. To analyze 
correlations in interference and diffraction, we need a quantum formalism of 
diffraction. In the early days of quantum electrodynamics(QED), it had been 
proved that the Maxwell equations which underpin diffraction remain true when 
the fields are quantizedl 15-17 ! . In quantum òptics the entire mode structure of 
the diffraction field is still determined by the Helmholtz part of the wave equa- 
tion. The role played by quantum mechanics is in determining the states of the 
diffraction modes from the states of the incident modes. However, no systematic 
approach in determining the states of the diffraction modes has been proposed. 
In this paper, we first solve this problem. By introducing the quantum Kirch- 
hoff boundary condition, we connect the states of the diffraction modes with 
the states of the incident modes by a diffraction factor. Then correlations of 
the diffraction modes with different kinds of incident optical fields are calcu- 
lated. Influence of correlations of the incident modes on the diffraction pattern 
is analyzed. The "ghost' diffraction is also explained based on this formalism. 

We consider Fraunhofer diffraction. This kind of diffraction is most im- 
portant. In Section 2, we introduce the equivalent scalar optical field and 
the quantum Kirchhoff boundary condition. The equivalent scalar optical field 
simplifies the problem of scalar diffraction, in which the variation of polariza- 
tion through diffraction is not considered. The quantum Kirchhoff boundary 
condition is equivalent in physics to the Kirchhoff boundary condition in clas- 
sical scalar diffraction yet overcomes the difficulty that the classical Kirchhoff 
boundary condition destroys the commutation relations of the field operators. 
In section 3, we obtain quantum formalism of Fraunhofer diffraction. The nor- 
mal characteristic functions of the diffraction modes are connected with those 
of the incident modes by a diffraction factor. From the characteristic functions, 
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correlation properties of the diffraction modes are analyzed in Sec. 4. In this 
section the diffraction pattern is also calculated with entangled incident states. 
An explanation of the "ghost" diffraction is proposed. 



2 The equivalent optical field and the quantum Kirch- 
hoff boundary condition 

In the diffraction problem the incident and diffraction optical fields are free. The 
free quantized electromagnètic field can be expanded into plane wave modes: 

E => > i\ — a-, e _> e V J+h.c, (1) 

where \x is polarization index and k ■ ~~e— ♦ =0. The annihilation and creation 

k fA, 

operators , at+ / ; satisfy the commutation relation 

k fe /lí 



^ A 4 fe /Li 



The frequency of the optical field remains unchanged through diffraction. 
So we only need consider fields with a definite frequency uj. That is, in the 
expansion (1) only the terms with k = ^ need be considered. Let k = 



k x , k y ,\J ^ — k'l — k^j and k has only two degrees of freedom k x ,k y (The 

symbols k below all have this meaning). The incident and diffraction optical 
fields are in a half space. Suppose the plane z = is the diffraction plane, then 
there may exist evanescent waves with a depression factor e~' fcz ' z at both sides 
of the diffraction plane. So the value domains of k x , k y are (— oo, +oo), i.e., k z 
can be imaginary. This is different from the plane wave expansions in the whole 
space. 

In scalar diffraction theory, the boundary condition at the diffraction plane 
is independent of the orientation of the optical field, and the variation of po- 
larization of the optical field through diffraction need not be considered. So 
we can introduce the following equivalent scalar optical field by neglecting the 
polarization index. 

where the box-normality of space has been used and S is the cross-section area 
of the box. The commutator (2) yields the following commutation relation of 
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the equivalent optical field at the diffraction plane z = 




(4) 



In scalar diffraction the equivalent scalar optical field can be in place of the 
real optical field. The diffraction problem is much simplified by introducing the 
equivalent scalar optical field. 

In classical scalar diffraction theory the Kirchhoff boundary condition states: 
the optical field remains unchanged through the diffraction aperture £ and de- 
cays to zero through the diffraction screenl 18 !. This boundary condition can 
not be used directly in the quantum case because the postulate that the optical 
field decays to zero through the diffraction screen destroys the commutation 
relations of the field operators. To keep consistent with quantum theory, we 
introduce the following quantum Kirchhoff boundary condition. The equiva- 
lent optical field e (x, y, z = 0) before diffraction is generally in a complicated 
entangled state and we use p(z = 0~) to represent its whole density operator. 
The quantum Kirchhoff boundary condition says: When passing the diffraction 
screen all modes of the field e (x, y, z = 0) (x, y € S — S) at the screen undergo 
such a strong dissipation that after the screen they are all in the vacuum state. 
At the same time, the modes of the field e (x, y, z = 0) (x, y G X) at the aper- 
ture undergo no dissipation at all. From quantum dissipation theory! 19 ' 20 ! , the 
total density operator p (z = + ) after diffraction is expressed as 



where the notation trss indicates trace of all modes at the screen. This 
boundary condition for scalar diffraction is equivalent in physics to the classical 
Kirchhoff boundary condition. Yet it is consistent with quantum mechanics 
as it results from the quantum dissipation theory. In next section we use this 
boundary condition to derive quantum formalism of diffraction. 

3 Quantum formalism of Fraunhofer diffraction 

In Fraunhofer diffraction the incident and diffraction optical fields are expanded 
into the plane wave modes and the role played by quantum mechanics is in de- 
termining the states of the diffraction modes from the states of the incident 
modes. Let a,—,, and represent the annihilation operators of the incident 



p( z = 0+) = trs-^p (z = 0") ® ]J \°)xy *y <°l . 



(5) 



(x,y)eS-X 



k 



k 
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mode k and the diffraction mode k , respectively. p ( a_>, ) is the density 



operator of the incident mode k and other incident modes are supposed in 
the vacuum state. First we derive the reduced normal characteristic function 
(b-^;^j of the diffraction mode b-g. Using Eq. (5) and the inverse trans- 



formation of Eq.(3) 



we get 



(6) 



X 



e k ■ e k 



= Tr{tr s _zp(z = 0-)® U |0>„ <0| 



(x,j/)es-s 



• cxp 



iÇ*-j= f s dxdye + (x,y,0) e l ( kxX+k y y ï • exp iÇ-^ J s dxdye (x, y, 0) e l ( k ^ x + k yv) j 



Tr { 



f> \ a t' 0/ 

k k q 



n |0)^, r# <0| exp l i? £ a+ i J s dxdye 



k' 



exp ^ i£ £ a-r*, | Je dxdye 



k' 



(7) 

where S and 5 represent area of the diffraction aperture and the whole diffrac- 
tion plane, respectively, and the notation Tr indicates trace of all modes. We 
define the energy transmissivity A as A = Its physical meaning is the ratio 
of the energy of the diffraction optical field to the energy of the incident optical 
field. The Fraunhofer diffraction factor / ( k ) is defined as 



/ k 



e -i{k X X + k yV ) dxdy _ 



(8) 



/ ( k lis normalized by 



e r (y) f (y) = i. 



(9) 



Eq.(7) is therefore simplified to 



X 



n) (^;ç)=x (n) [ a ^^f y- 



(10) 
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Eq.(lO) connects the reduced normal characteristic function of the diffraction 
mode b-g with that of the incident mode a-g, by a simple diffraction factor. 

Similar to the derivation of Eq.(lO), the total normal characteristic function 
of all diffraction modes j^} has the form 



#({^}>{^}) 



X 



(n) 



(11) 



The above results are obtained with the supposition that only the incident 
mode k ' is not in the vacuum state. If all the incident modes are in an 
entangled state, and we use Xt^ ({ a "^' } ' {^"^ ' }) ^° indicate its whole normal 
characteristic function. Eq. (11) can thus be generalized to 



n) 



~k ) 



(12) 



Eq. (12) determines the states of all diffraction modes from the states of the 
incident modes. It is a fundamental equation in the quantum formalism of 
Fraunhofer diffraction. 

The final result (12) is similar to the quantum description of the beam 
splitter. For the beam splitter, the input and output modes are linked by a 
canonical transformationt 21 ] 



6 2 



r 

-t 



ai 

«2 



(13) 



where a\,a2 are input operators and 61,62 are output operators. The parame- 
ters r and t should satisfy r 2 + í 2 = 1. From Eq. (13), we obtain the relation 
of the normal characteristic function between the input and output modes 



X in) (&i,&2;£i,&) = X (n> (ai,02;r£i -í&.ÍÉi + r&) 



» 



(14) 



Eqs. (14) and (12) are very alike in the form. However, some important differ- 
ences lie in their derivation. In diffraction the input and output modes cannot 
be put in a canonical transformation, which may be seen from the relation 



Y, V\f* [k - ~k) y/\f [k - ~k j = A < 1 . 



(15) 
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Only when the energy transmissivity A = ^ = 1 , i.e., when there is no diffrac- 
tion screen, the input and output modes can be linked by a trivial canonical 
transformation. So unlike Eq. (14), Eq. (12) is not a direct result of the input- 
output theoryt 22 ' 19 !. In the derivation of Eq. (12), the quantum Kirchhoff 
boundary condition plays an essential role. 

The general equation (12) can describe interference as well as diffraction. If 
there are two diffraction apertures Si, £2 , the diffraction factor / ^ k ^ simply 
becomes 



/ k 



e ~i(k x x+k y y) dxdy ^ 



(16) 



Si + S 2 JS1+S2 

where A = ^1^2. . When Si, S 2 tend to zero, Eq. (12) with this / ^ k ^ gives 
quantum description of the double-slit interference. 



4 Correlations in interference and diffraction 

4.1 correlations of the diffraction (or interference) modes 

In this subsection we consider correlations of the diffraction modes. Suppose 
all the incident modes except k ó are in the vacuum state. First we show that 
the diffraction modes are not correlated only when the incident mode k is 
in a coherent state. If the diffraction modes are independent, the decomposi- 

tion Xt ] (fe} ' fe}) =11 X [n) fe>^>) should hold - From E q- ( n ) this 



decomposition holds if and only if %( n ) 



has the following form 



X 



(n) 



i(€*a*+Ça) 



(17) 



i.e., the incident mode is in a coherent state. Under this condition, the diffrac- 
tion modes are not correlated and all in coherent states. With any other kinds 
of incident optical fields the diffraction modes are in an entangled state. 

The above discussion shows that the diffraction modes are generally corre- 
lated. In experiments correlation of the photon number is widely used, so we 
first calculate the correlation coemcient of the photon number of two diffraction 
modes .The correlation coemcient is defined by 



An_, An^ 

k 1 k 2 



Ara- 



2 



(18) 
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where n-+ (i = 1,2) denotes the number operator of the mode k j. After some 

h i 

calculation, from Eq. (11) we obtain 

F — 1 

V = r, (19) 

[(F n + h 1 -l)(F n + h 2 -l)}2 

where F n is the Fano factor of the incident mode k ' , i.e., 

K = ^7 — r~S ( 2 °) 



n- 



and hi in Eq. (19) is defined by 

1 



hi 



X 



f [ k i k 



(i = 1,2). (21) 



The relation between rj and F n is illustrated in Fig. 1 



Fig. 1 



If the incident mode is in a thermal state, F n = ( n_+, ) + 1 and r\ tends to 



fe 

its maximum value 1 with (n-gi ^ >> 1 • The correlation coemcient rj gets its 

minimum value — [(hi — 1) (J12 — 1)P 3 with the incident mode in a Fock state. 

Though 77 « 1 if the incident mode is in a thermal state with (n-gi ^ >> lj 
the diffraction modes are not correlated perfectly in this case. That can be seen 
from residual variance of the variables n_> and n_> in the linear regression. 

fel k 2 

The residual variance of the variable has the formi 23 ! 

fe i 



ri— , 

fe 



1? 



(22) 



°^-(F n + /n-l)(l-rç 2 ) 



where /?i and 02 are linear regression coefficients. Suppose h\ = h 2 and 
) >> 1, then 



Var ( - ftn^ - ft) « = 2 > . (23) 
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So in this case the residual variance is very large. In fact, the equation r\ — > 1 

under the condition ( n^, ) — > oo results from the infinite variance of n— > and 

\ fc o/ fe i 
n— > .We can not conclude from r? — > 1 that the diffraction modes are correlated 

perfectly. 

For the beam splitter, the correlation of the number operator of the output 
modes has the same form as Eq.(19). However, there are still some differences. 
First, the equation ^ = 1 holds for the beam splitter whereas in diffraction 
we have ^ < 1 . So for the beam splitter, the correlation coemcient of 
the output number operators can attain its minimum value -1 with the input 
mode in a Fock state. Second, in diffraction or interference correlations of many 
modes can be generated whereas the beam splitter is only used to prepare two- 
mode entangled states. 

Correlation coefncients describe correlation properties of a pair of specialized 
operators. Several approaches to the description of quantum entanglement have 
been proposed. In particular, Schlienz and Mahler interpreted the difference 
between the entangled state and the product state as the entanglement I 24 '. 
Suppose p is the density operator of the whole system and p a = tr^p , p^ = tr a p 
, where the subscripts a and b represent two subsystems. The Schlienz-Mahler 
measure is defined byt 24 l 



N(p) r 9i 
7 = V N(p) -l tr i ÍP ~ Pa ® Pb) \ ' (24) 

where N (p) indicates the dimension of the density operator p and 7 defined 
above satisfies < 7 < 1. However, the more recent papers distinguish quantum 
entanglement from classical correlations ' 25 ~ 30 ' . The entanglement is interpreted 
as the degree of inseparability. The entangled state is said to be inseparable if it 
can not be expressed as a mixture of product states of two subsystems. In this 
interpretation, the Schlienz-Mahler quantity 7 measures the total correlations 
rather than pure quantum entanglement. It is now believed that pure quantum 
entanglement can not be fully described by a single quantityt 29 ! . Bennett et. al. 
defined two quantities:! 26 ' 2 ^ "entanglement of formation" defined as the least 
number of shared singlets asymptotically required to prepare p by local oper- 
ations and classical communication, and "distillable entanglement" defined as 
the greatest number of pure singlets that can asymptotically be prepared from 
p by local operations and classical communication. And recently, Vedral et. al. 
introduced a new measure of entanglement t 30 !, which intèrprets the entangle- 
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ment as the minimum distance to all separable states. These measures have 
the desirable feature that their expectations can not be increased by local op- 
erations, but the disadvantage of being hard to evaluate because of the implied 
optimization. The question is still open in this direction. 

Though the Schlienz-Mahler quantity 7 in fact measures the total corre- 
lations, it is superior to the correlation coefficients, since it is not limited to 
specialized observables. In the following we use the Schlienz-Mahler quantity 
to analyze correlation properties of the diffraction modes. Before doing this, 
we first introduce the following lemma. 

Lemma. Suppose pi,p2 are two density operators of boson fields, and 
Xi^ (£) j (0 are normal characteristic functions of p\ and p 2 , respectively, 
then we have 

tr( Pl p 2 )= /^(fíx^-f)^ 1 '- (25) 

J TT 

Proof. If generalized functions (such as derivatives of delta functions) are 
permitted, the existence proof of P-functions of Boson fields has been given by 
Klauder and Sudarshant 31,19 ] . So tr (P1P2) can be expressed as 

tr (P1P2) = J P\ (cu) (oc\ P2 \a) d 2 ce 

(26) 

= irfPí (a) Q 2 (a) d 2 a, 
where P\ (a) and Q 2 (a) are P,Q-functions of the density operators p\ t p 2 re- 
spectively. The P,Q-functions are Fourier transformations of the normal and 
anti-normal characteristic functions. So Eq. (26) can be rewritten as 

tr(p lP2 )= Í^X { ; i) (0xí a) (-0, (27) 
J vr 

where (£) indicates the anti-normal characteristic function. Eq. (27) is 
equivalent to Eq. (25) . This completes the proof. 

We calculate the Schlienz-Mahler quantity 7 with a thermal incident optical 
field. From Eq. (11) the normal characteristic function of the diffraction modes 
k 1 and k 2 has the form 

(28) 

where (N) is the mean photon number of the incident mode. For thermal states, 
the dimension of the density operator N (p) — > 00. Eq. (24) together with Eq. 
(25) yields 

7 2 = X -T-2 —2 + —1 (29) 

x\X2 - 4:y z x\X2 - y z X1X2 
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where 



xí = 2 (N) A 
y = (N)X 



f( k 



+ 1 (i = 1,2), 



(30) 
(31) 

From Eq. (29) it is obvious that 7 tends to zero if (N) — > 00 or (N) — > 0. 



If the diffraction factor satisfies / ^ & ] 
27 = X2 = 2y + 1 and 7 is simplified to 

1 2 

7 



A: 



k 



+ 



Ay + 1 (3y + l)(y + l) ' (2y + l) 2 
The relation between 7 and y is illustrated in Fig. 2 



we have 



(32) 



Fig. 2 



From the figure we see the Schlienz-Mahler quantity 7 attains the maximum 
when y m 1.1. The maximum value is 0.25. With a larger or smaller mean 
photon number, the correlation of the diffraction modes decreases. 

4.2 Influence of correlations of the incident modes on the diffrac- 
tion (or interference) pattern 

To show the influence of correlations on the diffraction pattern, we consider 
the circumstance with two incident modes being in an entangled state. The 
entangled state is prepared by a beam splitter with the input mode in a Fock 
state. From Eq. (14) the normal characteristic function of the two incident 
modes has the form 



(n) 



1-7+/ ) 1-7+/ ) £-7+/ ) Ç-7+/ 



*1 



\n) , (33) 



where a denotes the input mode of the beam splitter. The diffraction pattern is 
shown by the n 
Eq. (12) gives 



shown by the mean photon number distribution ( ) of the diffraction modes. 



n- 



A / k 



k 



n- 



+ 



/ k 



k 



n- 



r[k -kAj[k 




+ h.c. 



(34) 



Ara 



rf [ k — k 1 ) — tf [ k 
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If the two incident modes are not correlated, i.e., if they are represented by the 
density operator tr~^, p <g) tr-g, p, the mean photon number distribution of the 



diffraction modes becomes 



n- 



rf k 



k 



+ 



tfik-k 



(35) 



So the two incident modes are superposed coherently if they are correlated and 
incoherently if they are not. 

The influence of correlations on the diffraction (or interference) is dramat- 
ically illustrated by the "ghost" diffraction (or interference) effect. In the ob- 
servation experiment of the "ghost" diffraction! 14 ], a light beam, which is gen- 
erated from spontaneous parametric down-conversion (SPDC) and consists of 
two orthogonal polarization components (usually called signal and idler), is split 
by a polarization beam splitter into two beams, and detected by two distinct 
pointlike photoncounting detectors for coincidence. A Young's double-slit or 
single-slit aperture is inserted into the signal beam. Surprisingly, an interfer- 
ence or diffraction pattern is observed in the coincidence counts by scanning the 
detector in the idler beam. Here we give an exact explanation of the "ghost" 
diffraction. For the SPDC, the output light is in a superposition of the vacuum 
and two-photon statesl 12 ! 



m = \0) + FY,a+ c+ |0> 

k k 

T 



(36) 



where the operators aí +/ and cL, i represent the signal and idler modes, re- 
spectively. The normal characteristic function of the signal and idler modes 

is indicated by x? ({ a f } > { c f } ! } ' }) ■ Then the si S nal H g ht 

meets a diffraction screen and the idler light remains unchanged. The second- 
order correlation coefficient between a fixed diffraction mode and arbitrary idler 
modes is to be measured. From Eq. (12), the normal characteristic function of 
the diffraction and idler modes has the form 

k - k ' ' ' 



(n) 



v l k 




With Eqs. (36) and (37), we obtain the second-order correlation coefficient 
between a fixed diffraction mode k and an arbitrary idler mode k ' 



(b. 



k' c ~k' 



6+ ò_>ci ; c^, 

k k k k 



X\F\ 



fi k - k 



(38) 
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It is directly proportional to square of the diffraction factor. The diffraction 
pattern occurs by fixing the diffraction mode and scanning the idler modes. 
Therefore, Eq. (38) explains the observation in the "ghost" diffraction. 
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Caption 1: The relation between the correlation coeficient 77 of the diffrac- 
tion modes and the Fano factor F n of the incident mode. We choose h\ = /12 = 
3. 

Caption 2: The relation between the Schlienz-Mahler quantity 7 and the 
mean photon number of the incident mode. y is expressed by Eq.(31). 
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